arXiv:math/0112143v4 [math.PR] 21 Jul 2003 


Growth fluctuations in a class of deposition 

models 


Marton Balazs 

(Institute of Mathematics, Technical University Budapest) 


Abstract 

We compute the growth fluctuations in equilibrium of a wide class of 
deposition models. These models also serve as general frame to several 
nearest-neighbor particle jump processes, e.g. the simple exclusion or the 
zero range process, where our result turns to current fluctuations of the 
particles. We use martingale technique and coupling methods to show 
that, rescaled by time, the variance of the growth as seen by a deter¬ 
ministic moving observer has the form \V — C\ ■ D, where V and C is 
the speed of the observer and the second class particle, respectively, and 
D is a constant connected to the equilibrium distribution of the model. 
Our main result is a generalization of Ferrari and Fontes’ result for simple 
exclusion process. Law of large numbers and central limit theorem are 
also proven. We need some properties of the motion of the second class 
particle, which are known for simple exclusion and are partly known for 
zero range processes, and which are proven here for a type of deposition 
models and also for a type of zero range processes. 

Resume 

On compute les fluctuations du grandissement dans I’etat d’equilibre 
d’une classe vaste des processus de decharge. Ces processus forment aussi 
bien un cadre pour quelques modeles des bonds voisins des particules, p. 
e. le modele simple exclusion on zero range, ou notre resultats deviennent 
des resultats sur les fluctuations du flux des particules. On utilise de 
methode martingale et des techniques des couplages pour presenter que le 
variance du grandissement, regradue par le temps et vu par un observateur 
qui avance determinement a une vitesse y, a la forme \V — C\ ■ D, oh C 
est la velocite de la particule de deuxieme classe, et D est une cons- 
tante connectee a I’etat d’equilibre du modele. Notre resultat principal 
est une generalisation du resultat de Ferrari et Fontes pour le modele 
simple exclusion. La loi des grandes nombres et le theoreme de la limite 
centrale sont aussi demontres. Nous avons besoin de quelques proprietes 
du mouvement du particule de deuxieme classe, qui sont connues pour 
simple exclusion et partiellement pour le modele zero range, et qui sont 
demontrees ici pour un type des processus de decharge et pour un type 
des modeles zero range aussi. 

Keywords: Current fluctuations; second class particle; coupling methods. 
MSC: 60K35, 82C41. 
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1. Introduction 


Stochastic deposition models can be used to obtain microscopic description of 
domain growths, e.g. a colony of cells or an infected area of plants. The fluc¬ 
tuation of the growth is itself of great interest. Moreover, these models are 
in close connection to interacting particle systems, where the particle diffusion 
corresponds to rescaled surface fluctuation. As it is shown below, an additional 
feature of deposition models is the possibility of handling antiparticles as well as 
particles in the particle representation of the process. It has been known [1] for 
the simple exclusion process, that the current fluctuation is in close connection 
to the motion of the so-called second class particle, and, divided by time, its 
variance vanishes for an observer moving with the speed of this particle. In this 
latter case, Prahofer and Spohn [2] suggest this quantity to be in the order of 
t2/3. 

In the present note we consider a wide class of one-dimensional deposition 
models, parameterized by rate functions describing a column’s growth depend¬ 
ing on the neighboring columns’ relative heights. By monotonicity properties of 
the rate functions, our models are attractive. For a treatment of these models 
in a hydrodynamical context, without using attractivity, see Toth and Valko 
[3]; Toth and Werner [4]. Following Rezakhanlou [5], we first show some con¬ 
ditions on the model in order to have product measures as stationary ones for 
the process. (By stationarity, we mean time-invariance in this paper.) Our de¬ 
scription is general enough to include the asymmetric simple exclusion process, 
some types of the zero range process, and a family of deposition models, which 
we call bricklayers’ models. In this general frame, we compute the growth fluc¬ 
tuations in order 0{t), hence generalize the result of Ferrari and Fontes [1]. In 
the computations we couple two processes, which only differ at one site. This 
is the position of the so-called defect tracer, or also called second class particle. 
We need law of large numbers and a second moment condition for the position 
of this extra particle. These have been established for simple exclusion [6], but, 
as far as we know, only L^-convergence is known for most kinds of zero range 
processes [5]. We prove ^"-convergence with any n for the defect tracer of the 
totally asymmetric zero range process and for our new bricklayers’ models via 
various coupling techniques. 

1.1 The model 

The class of models described here is a generalization of the so-called misan¬ 
thrope process. For —oo < w™'" < 0 and 1 < < oo (possibly infinite 

valued) integers, we define 

/ : = {z e Z : w™ - I < z < -f l} 
and the phase space 


n — {oj — (oJi)i£Z ■ COi £ 1 } — I . 

For each pair of neighboring sites i and i -|- 1 of Z, we can imagine a column 
built of bricks, above the edge {i, z -I- 1). The height of this column is denoted 
by hi- If w(t) G n for a fixed time t G K then uJi{t) = /ii_i(t) — hi{t) G / is the 
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negative discrete gradient of the height of the “wall”. The growth of a column 
is described by jump processes. A brick can be added: 


{iVi, oJi+i) —> {oJi - 1, LOi+i + 1) with rate r{uJi, oji+i). 

Conditionally on w(t), these moves are independent. See fig. 1 for some possible 
instantaneous changes. For small e, the conditional expectation of the growth of 
the column between i and z + 1 in the time interval [t, t + e] is Wi+i(t)) • 

e + o{e). 

The rates must satisfy 

r(u;™, •) = r-(-,a;““) = 0 

whenever either w™'" or is finite. We assume r to be non-zero in all other 
cases. We want the dynamics to smoothen our interface, that is why we assume 
monotonicity in the following way: 

(1) r{z + l,y)>r{z,y), r{y, z + 1) < r{y, z) 

for y, z, z + 1 G I. This means that the higher neighbors a column has, the 
faster it grows. Our model is hence attractive. 

We are going to use product property of the model’s stationary measure. 
For this reason, similarly to Rezakhanlou [5], we assume that for any x, y, z G I 


(2) r{x, y) + r(y, z) + r{z, x) = r{x, z) + r{z, y) + r{y, x), 

and for w™'" < x, y, z < -|- 1 

(3) r{x, y-1)- r{y, z - 1) ■ r{z, a; - 1) = r{x, z - 1) ■ r{z, y - 1) ■ r{y, x - 1). 

These two conditions imply product structure of the stationary measure, see 
section 1.3. Equation (3) is equivalent to the condition r{y, z) = s{y, z + l)-f{y) 
for some function / and a symmetric function s. 

At time t, the interface mentioned above is described by w(t). Let (/? : —> R 
be a finite cylinder function i.e. (p depends on a finite number of values of uii. 
The growth of this interface is a Markov process, with the formal infinitesimal 
generator L: 

(4) {Lip){u) = ^ r{uj^, uji+i) ■ [ip {..., Wj - 1, uji+i -1-1, ...) - V3(w)] • 

When constructing the process rigorously, problems may arise due to the 
unbounded growth rates. The system being one-component and attractive, 
we assume that, with appropriate growth condhions on the rates, existence 
of dynamics on a set of tempered configurations (i.e. configurations obeying 
some restrictive growth conditions) can be established by applying methods 
initiated by Liggett and Andjel [7] [8]. Technically we assume that fl is of full 
measure w.r.t. the canonical Gibbs measures defined in section 1.3. In fact this 
has been proved for some kinds of these models, see below. We do not deal with 
questions of existence of dynamics in the present paper. 
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1.2 Examples 

There are three essentially different cases of these models, all of them are of 
nearest neighbor type. 

1. Generalized exclusion processes are described by our models in case 
both and are finite. 

• The totally asymmetric simple exclusion process (SE) intro¬ 
duced by F. Spitzer [9] is described this way by = 0, = 1, 

r{uji, Wi+i) = Wi • (1 - uji+i). 

Here uji is the occupation number for the site f, and r(a>i, w^+i) is 
the rate for a particle to jump from site i to i -I- 1. Conditions (1), 
(2) and (3) for these rates are satisfied. 

• A particle-antiparticle exclusion process is also shown to de¬ 

monstrate the generality of the frame described above. Let = 
— 1, = 1. Fix c (creation), a (annihilation) positive rates with 

c < a/2. Put 

r(0, 0) = c, r(0, -1) = ^, r(l, 0) = |, r(l,-1) = o, 

and all other rates are zero. If Ui is the number of particles at site 
i, with Wi = —1 meaning the presence of an antiparticle, then this 
model describes a totally asymmetric exclusion process of particles 
and antiparticles with annihilation and particle-antiparticle pair cre¬ 
ation. These rates also satisfy our conditions. 

Other generalizations are possible allowing a bounded number of particles 
(or antiparticles) to jump to the same site. By the bounded jump rates 
and by nearest-neighbor type of interaction, the construction of dynamics 
of these processes is well understood, see e.g. Liggett [10]. 

2. Generalized misanthrope processes are obtained by choosing > 

—oo, = oo. 

• The zero range process (ZR) is included by w™™ = 0, = oo, 

Az-, y) = f(z) 

with an arbitrary / : Z"*" —> 1R+ nondecreasing function and /(O) = 
0. Here uJi represents the number of particles at site i. These rates 
trivially satisfy conditions (1), (2), (3). The dynamics of this process 
is constructed by Andjel [8] under the condition that the rate function 
/ obeys the growth condition \f(z + l) — f(z)\ < K for some K > 0 
and all z > 0. 

3. General deposition processes are the type of these models where 

u^min = —oo and = oo. In this case, the height difference between 

columns next to each other can be arbitrary in Z. Hence the presence of 
antiparticles can not be avoided when trying to give a particle represen¬ 
tation of the process. 
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• Bricklayers’ models (BL). Let 

r{z, y) ■■= f{z) + f{-y) 

with the property 

f(z)-fi-z+l) = l 

for the nondecreasing function / and for any z € Z. This process can 
be represented by bricklayers standing at each site i, laying a brick 
on the column on their left with rate f{—u)i) and laying a brick to 
their right with rate /(tOi). This interpretation gives reason to call 
these models bricklayers’ model. Conditions (1), (2) and (3) hold for 
r. Similarly to the ZR process, this model is constructed by Booth 
and Quant [11] only in case \f(z + 1) — f{z)\ is bounded in Z. 


1.3 Translation invariant stationary prodnct measures 

We are interested in translation invariant stationary measures for these pro¬ 
cesses, i.e. canonical Gibbs-measures. We construct such measures similarly to 
Rezakhanlou [5] of the following form. Fix /(I) > 0 and define 


( 5 ) 




for -I- 1. Then / is a nondecreasing strictly positive function. 

For I B z > 0 we define 


f{zy. := f{y), 

y^l 


while for / 9 2 < 0 let 


/W!: = 


1 


n fiy) 

y=z+l 

finally /(O)! : = 1. Then we have 

f{z)l-f{z + l) = f{z + l)l 

for all z € /. Let 

log fliminf (/(z)!)^^^i = lim log(/(z)) , if = oo 

\ Z—yoc / z^oo 


else 


and 


log (limsup (/(—) = lim log(/(—z)) , if w™'" = — oo 

\ z^oo J 2:^00 


, else. 


By monotonicity of /, we have 0 > 0. We assume 0 > 
parameter 0 G (0, 0 ), we define 


Z(»); = y 

ZGI 


e 


9z 


fizy- 


With a generic real 
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Let the product-measure have marginals 


1 e^^ 

( 6 ) ^ieiz)=ix^W: u;, = z}:=j^-j^. 

By definition it has the property 

^J■0iz + 1) _ e^ 
fig{z) f{z + l) 

which implies 


r{z + 1, 2 /- 1) 


fxeiz + l)fi0{y- 1 ) 

^J-9{z) fieiy) 


= r{y, z) 


due to (5) and (3). Hence stationarity of follows via (2). 

As can be verified, the expectation value q{9) : = Ee(a;i) is a strictly increas¬ 
ing function of 9. We introduce its inverse 9{q) and the function 


( 8 ) 


nig) : = Ee(g) {r(wi, w^+i)} , 


playing an important role in hydrodynamical considerations. For the SE model, 
the construction leads to the well-known Bernoulli product-measure with mar¬ 
ginals 

2x(l) = y{ui : LOi = l} : = Q, 
y(0) = /i{w : uJi = 0} := 1 — g 

with a real number g between zero and one (the density of the particles). In 
our notations, —g describes the average slope of the interface. 

For the particle-antiparticle exclusion process, the relative probability of 
having a particle or an antiparticle as a function of the rates goes as i/c/a, 
independently for the sites. The density of particles relative to antiparticles can 
be set by an arbitrary parameter. 

Both for the ZR process and for BL models, it turns out that / defined in 
(5) and / in the definition of the rates agree. 

It is not hard to show ergodicity of these models, which also implies ex- 
tremality of the invariant measures 

Proposition 1.1. The processes given in subsection 1.1, distributed according 
to their stationary measures (6), are ergodic. 

Proof. We need to show that any (time-) stationary bounded measurable func¬ 
tion defined on the trajectories of the process is constant a.s. By proposi¬ 
tion V.2.4 of Neveu [12], this follows once we see that any bounded function 
(p on H satisfying Pip = p is constant for /i-almost all w with the Markov- 
transition operator P. Hence ergodicity of the process follows if Lp = 0 implies 
p{ui) = constant for almost all w G We compute the Dirichlet-form 


- Tlg{p ■ Lp) = 



• b(- 


OJi — 1, CVi+l -|- 1, 


...) - pip)] 
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By positivity of the rates, this shows that assuming L ip = 0 results in 

ip{...,UJi- 1, Ui+I + 1, ...) = (p{uj) 

for almost all w G Consecutive use of this equation shows that^any function 
obeying Lip = 0 does a.s. not depend on any finite cylinder set in Q. Especially, 
for £ > 0 and a constant K G Ran((/?), the event 

{ip{uj_) G {K, K + s]} 

does not depend on any finite cylinder set. Hence by Kolmogorov’s 0-1 law, 
the probability of these events is zero or one w.r.t. the product measure p. 
Partitioning the bounded image of ip, this shows that this function is constant 
for almost all w. □ 


1.4 Results 


We start our model in a canonical Gibbs-distribution, with parameter 9. For a 
fixed speed value K > 0 we define 

:=hivtlit)-ho{0), 

the height of column at site \yt\ at time t, relative to the initial height of the 
column at the origin. For K < 0, we introduce 

:= h^vt]{t) - hoiO), 

which is the mirror-symmetric form of defined above for positive K’s. For 
V = 0 we write 

J{t) = J^°\t):=hoit)-hoiO). 

In particle notations of the models, is the current, i.e. the algebraic 

number of particles jumping through the moving window positioned at Vt, in 
the time interval [0, t]. We prove law of large numbers for this quantity: 

J{V) 

(9) lim -=E(r) —KE(a;) a.s. 

i^oo t 

We need law of large numbers and a second-moment condition for the po¬ 
sition Q{t) of the defect tracer (also called second class particle, see section 
3 for its definition) if one of the coupled models is started from its canonical 
Gibbs-measure: 

Condition 1.2. With initial distribution of u, weak law of large numbers 


( 10 ) 


lim Pg 

t—>oo 


Q{t) 

t 


c{0) 



= 0 


for a speed value C(9) and for any 5 > 0 holds, and the bound 


( 11 ) 



< K < oo 


is satisfied for all large t for the position Q{t) of the defect tracer. 
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Inequality (11) is obvious in case of bounded rates, since in this situation, the 
process |(5(t)| is bounded by some Poisson-process. 

Theorem 1.3 (Main). Assume condition 1.2. Then 

(12) lim ^ ^ Var,(u;o) = : Dj{e) 

t^OO t 

for any G K, where Varg stands for the variance w.r.t. fig. 

Theorem 1.4 (Central limit theorem). Assuming condition 1.2, 


lim Pg 

t—*oo 


v'DJm.Vi 



— oo 


i.e. J^^\t)!\/t converges in distribution to N{0, Dj{6)), a centered normal 
random variable with variance Dj{9) of (12). Tilde means here that the mean 
value of is subtracted. 

For the SE model, (10) is proven in [6]. It is shown there that 
lim = 1 — 2g a.s. 

t^OO t 

Condition 1.2 is satisfied by this law, hence theorem 1.3 gives 
lim g (1 _ g) 1(1 _ 2 p) - y|, 

and the central limit theorem 1.4 also holds. These results have been known for 
SE by Ferrari and Fontes [1]. 

For the ZR and BL models, we need a condition on the growth rates: 

Condition 1.5. For ZR and BL processes defined above, the rate function f is 
convex. 


For the ZR process, under this condition and assuming either strict convexity 
or concavity of 7i(p) defined in (8), more than (10), namely, -convergence is 
established by Rezakhanlou [5] with speed 


(13) 


C(0) = 


Varg (a 


As far as we know, the second-moment condition (11) has not yet been proven 
for this model. 


Theorem 1.6. For ZR and BL models satisfying condition 1.5 with initial 
distribution of gj, and for any n € , 


Q(t) 

t 


C(0) 


in L" 


where C{9) is defined in (13) for the ZR process, and 


(14) 


C{9) : 


2 sinh(0) 
Vare(a;) 


for the BL model. 



Hence under condition 1.5, condition 1.2 and thus theorem 1.3 and 1.4 hold 
for both ZR and BL models with C{6) defined in (13) and (14), respectively. 
As we expect by mirror symmetric properties of the BL model, the speed C(0) 
of the defect tracer is zero in case 0 = 0 in this model. 

Our methods do not rely on hydrodynamic limits. C{0) is a nondecreasing 
function for the totally asymmetric ZR process and BL model under condition 

1.5, see remark 5.10. This shows (non strict) convexity of the function 'H{g) of 
(8) for these models, since 

C{9{g)) = 

dg) 

after some computations, and 9{g) is also a monotone function. 

Proposition 1.7. Under condition 1.5, the function 7i{g) is strictly convex for 
the BL model. For the ZR process satisfying 1.5, linearity ofH{g) is equivalent 
to linearity of the rate function f on 1, which is the case of independent random 
walk of the particles. If this is not the case, then H{g) is strictly convex. 

This is an important observation for [13], since this property is only proved 
for small 9 values there. It is also remarkable for [5], where strict convexity is 
just assumed. 

We remark that rates for removal of the bricks can also be introduced to ob¬ 
tain a model with both growth and decrease of columns. In particle notations 
this represents possible left jumps of particles (or right jump of antiparticles, 
respectively). Therefore, not only the totally asymmetric case, but the general 
asymmetric case of particle processes (SE or ZR, for example) can also be in¬ 
cluded in the description. The extension of the proof of theorems 1.3 and 1.4 to 
this case is straightforward. However, the coupling arguments used to establish 
condition 1.2 for ZR and BL models in later sections are not applicable in case 
of brick-removal. 

We see that lim ^(*)) vanishes if we observe this quantity from the 

t^OO ^ 

moving position Vt = C{9)t, having the characteristic speed of the hydrodynam- 
ical equation. This has been known for the SE model with strongly restricted 
values of uji, and now it is proven for the class of more general models with pos¬ 
sibly uji G allZ also. The interesting question, of which the answer is strongly 
suggested for some models [2], is the correct exponent of t leading to nontrivial 
limit of Var(j(^^(t))/t^“ as t ^ oo. a is believed to be 1/3, in close connection 
to order fluctuations of the position Q{t) of the defect tracer. 

The structure of the paper is the following: after some definitions on the 
reversed chain, we begin with separating martingales from Var(J(t)) in section 
2.2. Then we proceed in section 2.3 by computing the generator’s inverse on the 
rates and then by transforming Var(J(t)) into nontrivial correlations. These 
correlations can be computed using monotonicity thus coupling possibilities of 
the model, this is done in section 3. This section also includes a technical 
lemma showing an interesting relation of space-time correlations to the motion 
of the defect tracer. After J{t), we deal with J^^'>{t), the growth in non-vertical 
directions in section 4. Our results are proven in this section, except for theorem 

1.6, which is proven in the last section for the totally asymmetric ZR process and 
for BL models. This last section includes the introduction of a new random walk 
depending on our processes, and new coupling techniques based on convexity 


9 



of the rate function /. As another consequence of these methods, this part is 
followed by a proof of strict convexity of the function 'H{q). 


2. The growth and correlations 

In this section we obtain a formula for Yar which contains only space- 
time correlations of iViitys as non-trivial expressions. 


2.1 The reversed chain 


The formal infinitesimal generator L* for the reversed chain is of the form 
{L*(p){u) = Wj+i) • [(p{.Ui + l, Wj+i - 1, ...) - V3(w)] 

on the finite cylinder functions. The rates r* of the reversed process w.r.t. 
can be determined by the equation 

Ee (V’(w) • Lip{^)) = Eg ((p(w) • L*y{^)). 

Proposition 2.1. For a;™" < z, y < 

(15) r*{z,y) =r{y, z). 

Note that the rates of the reversed process do not depend on the parameter 
9 of the original process’ distribution. 

Proof. Let ip, Lp be finite cylinder functions, and let X C Z be a finite discrete 
interval of which the size can be divided by three, and which contains the set 

{i ^ Ip, 01 {p depends on uji or on Wi_i} . 

Then the summation index i in the definition (4) of the generator can be run 
on the set 2. We begin by changing variables uJi, uJi+p. 


Ee (V’(w) • Lip{ui)) = 

= E,^{ r{uji, uji+i) ■ [ip{u)(p{.uj^-1, uji+i -I- 1, ...) - ipiaf) ■ V5(w)]} = 


iex 


iei 


= Ee r{uji + 1, Ui+i - 1) • 


pLe{ijji -I- l)pLe{uji+i - 1) 


yeicoi) y9{uii 


+ 1 ) 


X ip{..., LOi + I, Wi+i - 1, .. .)(/?(w)| - Ee I Wi+i)^ • • 

Since \I\ can be divided by three, we can apply (2) in order to show that 

'^r{uji, uji+i) = '^r{uj^+i, uJi) 
iei iei 

in the second term. By using (7) for the first term we finally obtain 


Ee (V’(w) • L(p{^)) = 

= Ee^ {r(wi+i, UJ^) ■ [ip{. ..,uji + l, Ui+i -1, ... )(p{u) - ip{u) ■ (/?(w)]} , 
iex 

which equals to Eg ((p(ui) ■ L*ip{ui)) by choosing r* according to (15). □ 
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Combining (7) with (15) leads to 


(16) 


r* (z, y) 


fX9{z + l)y.0{y- 1 ) 
Me (2) Me(j/) 


■ r {z + 1, y - 1), 


which is the natural formula suggested by considering conditional expectation 
values. 

In order to simplify notations, let 


r{t) :=r(wo(i), wi(t)), r*{t) : = r*(ivoit), wi(t)). 


2.2 Preparatory computations 

For a quantity A(w) with E|A| < 00 , let 

A = A(w) := A-EA. 

t 

Lemma 2.2. Var(J(t)) = t E(r) +2 f (t — v)E (r(v) r*(0)) dv. 

0 

Proof. By definition, E(J(t) | w(0)) = <r(0) + o(t), hence 


M{t) : = J{f) — J r{s) ds 
0 

is a martingale with M(0) = 0. Using this, 

(17) Var(J(t)) =EM(t)2 + 2E I^M(t) yf(s)dsj +E ^j¥{s)ds 
Due to E llil(O)) = tr{0) + o{t), the process 

t 

N{t) : = M{tY — J r{s) ds 
0 

is also a martingale with iV(0) = 0. Hence 

EM{tf = tE(r). 

Using the martingale property of M, the second term of (17) can be written as 

i i 

2 J E{M{t)r{s)) ^s = 2 J E(M(s)f(s)) ds. 

0 0 

Simply changing the limits of integration in the third term of (17), we have 


E 
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These calculations lead to 


(18) Var(J(t)) = t E(r) + 2 J E |^r(s) |^M(s) + J r('u) duj J ds = 

t 

= t E(r) + 2 y E ms) J{s)) ds. 


In order to handle E(r(s) J(s)), we introduce the quantity corre¬ 

sponding to J in the reversed model by 

: = J{s) — J{s — u) (s > u > 0). 

This is the number of bricks removed from the column in the reversed model 
started from time s. As in case of J{t), a reversed martingale can be separated 
by 

u 

- J r*{s-v) du. 

0 

For this reversed object, = 0 and E | = M(®)*(s — 

t) if 0 < s — t < u, where T stands for the natural filtration of the (forward) 
process. In view of this, 


E (f(s) J(s)) = E f(s)E(j(®)*(s) 




[s, oo) 


= E I r(s) J r*{s — v)dv \ = J E, (r(v) r*(0)) dv, 
Vo / 0 

where in the last step we used time-invariance of the measure. Using this result, 
we obtain 

t 

Var(J(t)) = t E(r) +2 J (t — v) E (r(u) r*(0)) du 


from (18) by changing the order of integration. 


□ 


2.3 Occurrence of space-time correlations 

In this subsection we denote r{uji, uii+i) and r(wi, uii+i) by and r^, respec¬ 
tively. For fc e Z, let 


dfc . U ^ I , 

be the k-th. coordinate of U. Then 


(19) 


(Ldfc) (w) = Tk-i - Tk and 
{L*dk){ui) = -r^_i +rl , 
where L* is the infinitesimal generator (2.1) for the reversed process. 
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Lemma 2.3. For 0 < a < 1 the expressions 

OO OO 

( 20 ) (fia ■ = '^a^~^dk ipa ■='^a’‘d-k 

k=l k=0 

exist a.s., and 

\\Ta{Lipa){ui) = — lim (L'!/)q,)(w) =r, 

Oi —»^1 Oi —»^1 

\Ym{L*'ipa){u£) = — lim (L*(/9Q,)(w)=r* 

Oi — Oi —»^1 

in L^. 

Proof. The a.s. existence of the sums above can be easily shown by using the 
Borel-Cantelli lemma for the sets 

A„ : = {^ : |w„| > n} . 

We show the first equation for ipa- By (19) 


(21) (T</?a) (w) = ro + (a - 1) ^ = 

k=l 

OO OO 

= ro - Er + (a - 1) ^(r^ - Er)a''“^ = ?o + (a - 1) y^?fcQ;^~b 
fe=i fe=i 

By independence of coi and ujj for i ^ j, E(n • rfc) = 0 if |Z — /cj > 1 and 
|E(r) • rfe)| < E(r) • n) = ||r|| 2 , if |fc — ^| = 0 or 1. Hence the L^-norm of the 
second term on the right-hand side of (21) tends to zero as a ^ 1: 

2 

<(a-l)2^||ffc||2a2fc-2 + 

2 k—l 

( -\\2 

+ 2(0 - 1)2 ^ WkWW^-^ = ^f^lFll2(l + 2a-^) -- 0. 

fe=l 

The proof of the other three equations is similar. □ 

Now we can compute the integrals in our expression for Var(J). 

Theorem 2.4. 

OO 

Var(J(t)) =tE(r) -2tE(r*(0) • wi(0)) + 2^nE(wo(0)w„(t)) = 

n—1 

OO 

(22) =tE(r) -h2tE(r*(0) • wo(0)) + 2^nE(wo(0)w-„(t)). 

n—1 

As can be seen in the next session, the sums on the right-hand side are 
convergent. 


{a - 1) '^rua^ ^ 

fc=i 
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Proof. Using convergence stated in lemma 2.3 and Cauchy’s inequality, we 
rewrite the integral in the result of lemma 2.2. We can write P(0) instead of 
r*(0) there, since 'Ei{AB) = 'E,{AB) if both sides exist. 


(t — u)E(r(u) r*(0)) du — lim f {t — v)E{LLpa{v)r* (O)) dv < 

a^l J 
0 
t 

< lim [ {t — v) a/E ([r(v) — Lipa{v)]'^) ■ E(r*(0)2) dv = 
a^l J 
0 

t 

= lim [ {t — v)dv = 0, 

a^l J 


hence we can apply integration by parts: 

t 

Var(J(t)) = t E(r) + 2 lim [ {t — v)E {Lipc,{v) r* (0)) du = 

0 

t 

(t — v) —E((^q(u) r* (0)) du = 
dv 

0 

i 

= t E(r) — 21 lim E((Pq( 0) P(0)) + 2 lim f E((^q,(v) P (0)) dv. 

0 

The last integral here can be transformed in the same way, using lemma 2.3 
again: 


z 

E((?„(u)P(0))du = j E(^„(0)P(-u))du = 

0 

t t 

= lim f E(ipaiO) L*^l;.y{—v)) dv = lim f ——E{ipaiO)ip-y(—v))dv -- 
1^1 J 7^1 J dv 


= lim E((Pq,(0) — lim E((Pq,(0) '^^-,,(0)). 

7.^1 7^1 
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Hence with definitions (20), the variance of J{t) can now be written as 


Var(J(t)) = tE(r) - 2 Him E(^„(0) P (0))+ 


a^l 


+ 2 lim E{ipa{0) — 2 lim E((^q,(0) ^(O)) = 

a:,7—Q:,7— 


= tE(r)-2t ImE Wfc(O) P(0)j + 

/ oo oo \ 

+ 2 Jim ^ E I ^ a'^~^ Wfc(O) ^ 7' uj-i{-t) j - 

/=o / 

/ OO 00 \ 

2 lim E I Q;^“^Wfc(0)y^ a*w_/(0)). 

Q:,7—\ / 






/=0 


Using product property of the measure at time t = 0 and the fact that r* 
depends only on ojq and oji, most of our expressions become simple (recall that 
all quantities with tilde are centered random variables): 


Var( J(t)) =tE(r) - 21 E(5i(0) P(0)) + 2 EE E(wfc(0) w_/(-t)) - 0 = 

k^l1^0 

00 

= tE(r) - 2tE(wi(0)P(0)) + 2yjnE(w„(t)wo(0)). 

n—1 

In the last step, we used translation- and time-invariance of the measure. 

We needed Lipa —> r and L*'ipa. P in so far. The properties —Lipa ^ 
and —L*ipa r* can be used in a similar way to prove the second equation of 
the theorem. However, we need both tpa and ipa- using only one of them would 
have lead to a divergent sum in the last step. □ 

The first two expressions of formula (22) can be computed easily. The diffi¬ 
culty is in determining the space-time correlations E(wo(0) Wfe(f))- In order to 
do this, we use coupling technique. 


3. Coupling and correlations 

In this section, we show how to couple a pair of our models, with the help of 
the so-called second class particles. We can use second particles to compute our 
expressions containing space-time correlations. 

3.1 The basic coupling 

We consider two realizations of a model, namely, C and rj. We show the basic 
coupling preserving 

(23) Q{t) > 77,(t), 

if this property holds initially for t = 0. We say that n = (i{t) — rii{t) > 0 is the 
number of second class particles present at site i at time t. During the evolution 
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of the processes, the total number of these particles is preserved, and each of 
them performs a nearest neighbor random walk. 

The height of the column of C (or rj) between sites i and t + 1 is denoted by 
gi (or hi, respectively). (These quantities are just used for easier understanding, 
they are not essential for the processes.) Let gi t (or hi '[) mean that the column 
of C (or the column of rj, respectively) between the sites i and z + 1 has grown 
by one brick. Then the coupling rules are shown in table 1. Each line of this 
table represents a possible move, with rate written in the first column. In the 
last column, (or ,r^) means that a second class particle has jumped from z 
to z + 1 (or from z + 1 to z, respectively). This coupling for the SE model is 
described (with particle notations) in Liggett [14], [10] and [15]. The rates of 
these steps are non-negative due to (23) and monotonicity (1) of r. These rules 
clearly preserve property (23), since the rate of any move which could destroy 
this condition becomes zero. Summing the rates corresponding to either gi ] or 
to hi t shows that each C and zy evolves according to its own rates. It would be 
possible to couple models possessing rates for removal of bricks as well. 


3.2 Correlations and the defect tracer 

We introduce the notation G 12, a configuration being one at site z and 
zero at all other sites. Let w be a model distributed according to and 
C(0) = a;(0) -I- 5 q, i.e. we have a single one second class particle between ^ and 
w, initially at site 0. In order to avoid confusions, we call this particle the defect 
tracer. According to the basic coupling, this single defect tracer is conserved 
for any time t: 

(24) C(^) = + ^Q(t) 

The quantity Q(t) is the position of the defect tracer, performing a nearest 
neighbor random walk on Z. 

In this subsection we consider the process (w(t), Q{t)), the model distributed 
according to the Gibbs measure g and the random walk Q{t) connected to it 
with (5(0) = 0. Using condition 1.2, we prove theorem 1.3 for U = 0. We begin 
with a technical lemma, showing how to make use of the defect tracer. 

Lemma 3.1. For the pair (w(<), Q{t)) defined above, and for a function F : 
/ ^ R with F(a;"*“^) = 0 and with finite expectation value ^ T'(z) g{z), 


(25) 


E 



FjcvoiO) - 1) gjiUoiO) - 1) 
yz(wo(0)) 


F{toom 


= E(l{g(<) = 


rz}F(a;o(0))). 


Proof. We take conditional expectation value of (24): 


(26) E (C„(t) I wo(0) = z) = E (w„(t) I wo(0) = z) -f P {Q{t) = n \ wo(0) = z). 


Initially, C(0) = w(0) -I- 6q. Therefore, ( itself is also a model with initial 
distribution g, except for the origin. Hence 

E (C„(t) I a;o(0) = z) = E (C„(t) | Co(0) = z + 1) = E (a;„(t) | a;o(0) = z + 1), 
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and (26) can be written as 

E (w„(t) I a;o(0) = z + 1) - E | wo(0) = z) = P {Q{t) = n \ wo(0) = z). 
We multiply both sides with F{z) /r(z) and then add up for all z € / to obtain 

E (w„(t) I a;o(0) = z) • (F(z - 1) ^l{z - 1) - F{z) ^(z)) = 

ZGI 

= X! P (Qit) = n I Wo(0) = z) • F(z) ^J,(z). 

zei 


Here we used that = 0 and we write — 1) = 0. We know that 

P(a;o(0) = z) = /r(z), hence the proof follows. □ 

Corollary 3.2. We use the convention that the empty sum equals zero. Let 

g{z) ■■ = z-'^yp,{y). 

yGl 


For n € 7 j, 


E(iio(0)S„(t)) = E(^l{g(t)=n}. ^ 

Proof. By the previous lemma, our goal is now to find the correct function F, 
for which 

F(z —l)^(z —1) 

--F(z) = g{z) = z - 2_^yg{y) 

is satisfied. By inverting the operation on the left side, we find 

, .max 

Liv) 


hjz) \ 


F{z):= 9{y) 

y=z+l 


/i(z)' 


This function satisfies the conditions of the lemma. Using (25), 

E(w„(t) wo(0)) = E(a;„(t) wo(0)) = E(w„(t) • 5(a;o(0))) = 
■F(u;o(0)-1)m(wo(0)-1) 


. 

= E(l{Q(t) =n}E(a;o(0))) = 

l{g(t) =n}- g{y) 

y=uioiO) + l 


Ky) 
m(wo(o)) j ■ 


□ 

Now it becomes clear that we need to know something about the motion 
of the defect tracer, f and w can not be started together from their original 
stationary distribution due to the initial difference between them, present at the 
origin. We could follow our defect tracer. Knowing a measure stationary as seen 
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from site Q{t) for all time t would help us to state the law of large numbers for 
the Q(t) process. In general, we don’t know such a stationary measure which 
has the same asymptotics far on the left and far on the right side. It is shown 
in [13], that under some weak assumptions for BL models, this measure can not 
be a product-distribution. (Instead, a shock-like stationary product-measure is 
described there for certain type of rates, under which the slope of the surface 
differs on the left side from that on the right side.) 

For SE and some types of ZR processes, law of large numbers (10) is known. 
This law and the second moment condition (11) for BL and ZR models pos¬ 
sessing convexity condition 1.5 are proven in section 5. As shown in the next 
theorem, this allows us to do further computations on the space-time correla¬ 
tions of the models. We need the following properties of the canonical measure: 

Lemma 3.3. (i) The sum 


H 9{y) 

z&I y=z+l 


Ky) 


is convergent, and 
(ii) the sum 

max 

UJ 

E E giy)Tiy) = Var(wo) 

zeiv=z+i 

is convergent and the equality holds. 

Proof. For 9 e (9, 9), the tails of the measure have exponential decay. 
Hence the convergence in both expressions holds. The identity in (ii) is straight¬ 
forward and is left to the reader. □ 


The next lemma shows the essential connection of the defect tracer to space- 
time correlations in the model. 

Lemma 3.4. Assume condition 1.2 with speed value C. Let B{t) he a real¬ 
valued function with limj^oo H(t) = R G R, ni, n 2 G Z, A S M, Vi < V 2 in 
R U {— 00 , 00 } and the real interval V : = [Vi, V 2 ]. If either 

(i) C fz Vi, V 2 , or 

(ii) C G R and A - C = —B 
holds, then 


LtV2j+n2 

^Ihn ^ -k H(t)^ • E(wo(0) w„(t)) = 

n=('tVil+ni 

= {AC+ B)- IjC G V} • Var(a;o), 

where Var(a;o) is the variance of ujq w.r.t. the canonical Gibbs-measure. 
Proof. We define V* by 


V' : = 


El 


t' 


V 2 + 


n2 
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By corollary 3.2, 


[tV2i+n2 

(27) ,1m Y. + •E(wo(0)S„(t)) 


n=[tVi'l+ni 
ltV2\+n2 


^max \ 

^ + E I l{Q(t) =n} • Y 9{y) I 


n=['tVi]+ni 


y=wo(0) + l 


( ^max 

+ B(t)y i{g(i)/i e v‘} • ^ 5(2/) 


m(^o(0)) / 

y-{y) 


y=wo(0) + l 


m(wo(0)) 


= ^lim Y ^(+ ^(^)) ■ 1{Q(^)A G V*} • l{u;o(0) = z}] x 
ze/ \ \ / / 


XI 5(2/) 


y=z+l 


A^(2/) 

H{z)' 


We show that the limit and the summation can be interchanged in this expres¬ 
sion. We use Cauchy’s inequality to obtain 


E M 


MW 


-b B{t) ■ i{Q{t)/t e V‘} • 1M(0) = z} 


< 


< 


\ 


E u 


MW 


B{t) 


( Q{t) 


€ and a;o(0) = z I < 


< K' 


f Q{t) 


G V‘ 


wo(0) = y • yYY < K' ■ a/M) 


for some constant K' by (11). Since g{y) is monotone in y and 


, .max 

to 

Y 5(Wm(2/)=0: 


the sum 

, .max 

to 

Y 5 ( 2 /) m(2/) 

y=z+l 

is non-negative for any z € I. Hence we can bound from above the absolute 
value of the terms in (27) for each 2; G / by 


K'- 


max 

to 


Y 9iy) 

y=z+l 


yjy) 

\/MW’ 


and the sum 


^max 


Y^'- Y 9{y) 

z&I V=Z+1 


m(2/) 
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is convergent by lemma 3.3. Using dominated convergence, we write 
[tV2]+n2 

(1™ ^ • E(wo(0) w„(t)) = 

n=\tVi'\+ni 

= ,1™ E (a^ + Bit)] ■ l{Qit)/t G V*} • l{u;o(0) = 4 ) x 

zGl \\ / J 


X 


We introduce the set V‘ : = V*n,8e(C'), where for e > 0, BdC) 
R. Hence V‘ = V* U (V* \ 8^(0)): 


9{y) 


yjy) 

^i{zy 


V=z+1 

= (C-s, C+e) C 


[tV2j+n2 

(28) ^lim ^ • E(wo(0) w„(t)) 

n=[tVi'l+ni 


= /iz E [ [h® + B(t)) ■ l{Q(t)/t G V4 • l{^o(0) = z}) X 
zGl \\ / J 


giy) 


y=z+l 


yjy) 

giz) 


Y , 1 ™ E + Bit)^ ■ l{Qit)/t G V‘ \ B,iC)} ■ l{a;o(0) = 4 ) x 


Y 9iy) 


y=z+l 


B{y) 

Kz)' 


(28) contains two terms. We use Cauchy’s inequality on the second term as we 
have done before: 


E ( ( + B{t) ) • l{Q(t)A G V‘ \ BeiC)} ■ l{u;o(0) = z} 


< 


<K' -Jp 




G V* \ BeiC) and a;o(0) = z < 


<K'-Jp 


Qit) 


t 


i ByC) ^ 0 


as f ^ 00 by the law of large numbers (10). Only the first term of (28) remained, 
for which we write 


(29) 


\tV2\+n2 

^Ihn Y + H(t)^ • E(wo(0) w„(t)) = 

n=rtVi]+ni 

= lim (A - C + B{t) + 0{s)) ■ P iQ{t)/t G Vl and wo(0) = z) x 

t^OO ^ ' 

Z^I 


X 


^max 


Y 3{y) 

y=z+l 


yjy) 

giz)' 
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We have three possibilities. 

(i) If C G V, C ^ Vi, V 2 , then for small e and large t, V* = Be{C), and by (10), 
lim P (Q(t)/t G Vt and wo(0) = z) = 

t^OO ^ ' 

= lim P G Be{C) and wo(0) = z) = P(wo(0) = z) = fJ,{z). 

t—^OO 

Hence we can continue (29) by 


ltV2j+n2 

X! (J^A +■E,{u;o{0)ujn{t)) = 

= ^ iiin{A-c + B{t) + Ois))- g{y)Ky)^ 

t^oo 

z&i y=z+i 

^max 

(A-C + B)- E E giy)yiy) = {A-c + b)- Var(wo) 

zg/ y=z+l 

as £ ^ 0. The last equality is a result of lemma 3.3. 

(ii) If A - C = —B, then the right-hand side of (29) tends to 0 (e) as t ^ 00 for 
all e > 0, hence is zero in this limit. Here we used that 

P {Q{t)/t G Vg and wo(0) = z) < P(wo(0) = z) = iJ,{z), 


and that 

^max 

E E 9iy)yiy) 

zGl y=z+l 

is convergent. 

(iii) In case C ^V, then V* is empty for £ small and t large enough, and hence 
the right-hand side of (29) is zero. 

The result of these three cases completes the proof the lemma. □ 

Now we are able to compute limj^oo for V = 0. The proof of the 

general formula (12) requires some more computations in the next subsection. 

Theorem 3.5. Assume condition 1.2 with speed C. Then 

lim = E(r) - 2E(r*(0) ■dA{0)) + 2C+ ■ Var( wq) = 

t^OC t 

(30) = E(r) -f 2 E(r*(0) • wo(0)) -b 2 C" • Var(a;o). 

Here 0 < is the positive or the negative part of C, respectively. 

Proof. We consider the result of theorem 2.4. Dividing (22) by t and taking the 
limit t ^ 00 allows us to use the result of lemma 3.4. For the hrst equality of 
(22), we use this lemma with parameters Vi = 0, V 2 = 00 , ni = 1, n 2 = 0, A = 
1, B{t) = 0. Then we obtain 

lim = E(r) - 2 E(r*(0) • wi(0)) -b 2 C • IjC > 0} • Var(wo)- 

t—^'OC t 
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For the second equality of (22), we rewrite the sum as 

li^ Va^^E(^) + 2E(r*(0 )-wo(0))-2 V nE(io(0) 

t^OO t ^' 

n— — oo 

in order to use lemma 3.4 with parameters Vi = — oo, V 2 = 0, ni = 0, n 2 = 
— 1, A = 1, B(t) = 0. Hence 

lini = E(r) + 2 E(r*(0) • wo(0)) - 2 C • IjC < 0} • Var(a;o), 

t^oc t 

which proves the second equality of the theorem. □ 

We obtained two formulas for the variance of J{t). If the characteristic speed 
C exists, then we can compute it by subtracting the two lines of (30). 

4. The growth in non-vertical directions 

We have examined so far Var(J(t)), the growth fluctuation of a fixed column, 
i.e. the fluctuation of vertical growth. In this section we deal with the growth 
fluctuation of the surface in equilibrium, but considered in a slanting direction, 
namely, Var(J*^'^)(t)). From now on, we assume without loss of generity /io(0) = 

0 . 

Proof of (9). By definition ujj{t) = hj-i{t) — hj{t), we have 

i 

(31) hi{t) = ho{t) 


for any site i > 0, hence for V >0, 

[Vt] [Vt] 

= hyvt\{t) = ho{t) - = ho{t) - 

j—1 L J 

By ergodicity, the first term has the limit Er a.s. when divided by t. The second 
term is \yt\ times the average of an increasing number of different iid. variables. 
These variables have finite moments, hence the fourth-moment argument (see 
e.g. [16, Theorem 7.1]) is applicable with the discretization series : = n/V to 
show that 



This shows (9) for the limit taken along the subsequence t„. For any t € 
there is a unique rit € for which t < tnt+i, and 

is the number of bricks laid on column rit in a time interval shorter than 1/E, 
hence dividing it by t leads a.s. to zero in the limit. Therefore (9) holds for the 
limit of as well. Similar computation works for E < 0, and finally, the 

case E = 0 is trivial. □ 


22 



Now we consider the fluctuations (with tilde meaning the mean value sub¬ 
tracted). 

(32) Var(j('^)(t)) = E | (j(^)(t) - Ej('^)(t))^ 

= E |(hLvtj(t)) I = E hLytj(<) - hLvtj(O) -l-hLvtj(O)^ | = 

= E I {hyvt\ {t) - hyvt\ (0)) "I - E I {hyvt\ (0)) "| + 2 E {hyvt\ {t)hyvt\ (0)) • 

By translation-invariance, the first term is Var( J(t)), computed in the previous 
sections. By (31) and by product structure of the measure, the second term of 
the right-hand side of (32) is 

-E|(hLyij(0))'| =•E(iio(0)2) =-LWJ • Var(u;). 

The limit of the third term divided by t in (32) is computed in the following 
two lemmas: 

Lemma 4.1. For E > 0, 

(33) E (hyvtiit)hyytiiO)) = 

ivt\-l -1 

= ^ ([ytj - n)E(w„(t) Wo(0)) + ^ nE(w„(t) wo(0)). 

n— — oo n— — oo 

Proof. Using (31) again, 

(34) E(hLwj(t)hLvtj(0)) = 

/ LV^iJ \ fivt] [Vt\ 

= -E hoit) ^ S,(0) + E ^ w,(t) ^ 5,(0) 

V 1=1 / 1=1 

LVij ivtj ivt\ 

= - ^ E(ho(t)5,(0)) + E E E(5.(t)5,(0)). 

1 = 1 i=l 1=1 

A martingale 

t 

H{t) : = ho{t) - j ro{s)ds 
0 

with H{0) = 0 can be separated in order to show that 

t 

E(ho(f)5,(0)) =E(il(t)5,(0)) + J E(ro(s)5,(0))ds = 

0 

t 

= j E(fo(s)5i(0))ds. 

0 
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Now we use an argument very similar to the proof of theorem 2.4. By lemma 
2.3, the L^-convergence 

- lim {L'ipa){u) = ro 

can be used to replace our integral: for j > 1 we continue by 

t t 

E(/io(t) S3j(0)) = J E(ro(s) S;j(0)) ds = — lini J 'E{Ljpa{s)ujj{0))ds = 

0 0 
t 

= - lim [ ^E(V'a(s)wj(0))ds = E(V’a(0^i(0)) - E(V'a(0)wj(0)). 
a—»-l J as 
0 

Using definition (20) of ipa and product structure of the canonical measure, 


E(/io(t)w,(0)) =E 


-E 


i=0 


y]w_i(o)wj(o) = 


i=0 


oo 0 

2=0 i— — oo 


Combining this expression with (34) leads to 


0 lvt\ [Vti ivt\ 

E (h^vtiit)hvti{0)) = E E E(i.(t)w,(0)) + E E E(S.(t)ii,(0)) = 

i— — oo j — 1 i—1 j — 1 

[Vt\ [Vt\ [Vt\ [Vt\ 

= E E= E 

i— — oo j—1 i— — oo j—1 

by translation-invariance. Changing the summation indices leads to the proof 
of the lemma. □ 

Lemma 4.2. Assume condition 1.2. Then for U > 0, 

yin t E (0)) = (E - C+)+ • Var(u;). 

Proof. We use lemma 3.4 for the two terms on the right-hand side of (33). For 
the first one we set Vi = — oo, V 2 = V, ni = 0, 712 = —1, A = —1, B{t) = 
[Vt\/t, while for the second term in (33) we put Vi = — 00 , V 2 = 0, ni = 
0, n 2 = —1, A = 1, B{t) = 0. One can easily check that for any C G K. and 
U > 0, one of the cases (i) or (ii) of lemma 3.4 apply. Consequently, we obtain 

y™ \ ^ {^Vvt\ {i)\vt\ (0)) = 

= [{V -C)-1{C<V} + C- 1{C < 0}] • Var(w) = (V - C+)+ ■ Var(w). 


□ 
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Now we divide equation (32) by t and take the limit t ^ oo. We use the 
result of lemma 4.2 to obtain 

(35) ta YEWm + [2,v _ - n . Var(.) 

i—>cxo t t—*oo t 

for > 0. 

For V < 0, we proceed as we did above with for positive F’s. The only 
important difference is using ipa instead of —ipa in the proof of lemma 4.1. The 
result of a similar lemma for 14 < 0 is 




= - ^ i\Vt] - n)E{uJn{t) uipjO)) - nE(a;„(t)cjo(0)). 

n=|'Vi]+l 


n—1 


Therefore, lemma 3.4 is applicable in a similar way as in lemma 4.2 above. The 
result of this application is 




= (E + C-)- • Var(w). 


Computing Var(J^^)) for C < 0 as we did in (32) leads then to 

(36) lim --^41 = lijjj -(—(41 _|_ [2(V + C ) + V] ■ Var(a;). 

t — ^OO t t^OO t 

Now, assuming condition 1.2, we can prove (12) by the result of theorem 

3.5. 

Proof of theorem 1.3. All time arguments of our variables for this proof are 
thought to be zero without mentioning it. By (16), 


E {r* ■ (io - wi)) = E (r* • (wq - ^i)) = 

IT ,1 + 1) - 1) ^ 

- E (^r(u;o + 1, cui - 1) ■ 

E(r-(wo-u;i))-2E(r) = E(r-(wo-wi))-2E(r) = -E(r*-(wo-wi))-2E(r), 
we used (15) in the last step. Hence we obtain 

E(r*-(wo-wi)) = -E(r). 

We have two formulas for the variance Var(J(<)) by theorem 3.5, which are 
used together with (35) and (36) to obtain 

lim ^ _ 2 (|y _c\+C)- Var(w) = 

i^oo t 

= E(r) + 2 E(r* ■ uq) + {\V - C\ - C) ■ Var(w). 

We take the average of these two formulas: 

lim - (111 = E(r) + E (r* • (wq — wi)) + |y — Cl • Var(u;) = 

t^oo t 

= lE-CI • Var(w). 


□ 
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Now it is easy to prove central limit theorem for 
Proof of theorem 1.4- We introduce the drifted form of by i G Z: 

jP(t) ■■ = hoti+z(t) - h^iO) 

for C > 0, and 

jF\t) ■=hict]+zit) - h^iO) 

for C < 0. Due to translation-invariance, the distribution of this quantity is 
independent of i. Hence by (12), for C > 0 and H > 0, the variance of 

(^) = ^ivt} (t) - hivt}-ict\ ( 0 ) = - hivt\-ict\ ( 0 ) 

is o(t) as t oo. Thus it follows that we only need central limit theorem for 
(0), which is, by (31) and by /io(0) = 0, the sum of \ \yt\ — \ Ct\\ 
number of iid. Wi(0) variables with finite moments. Hence the theorem follows 
for H > 0, C> 0. For H > 0, C < 0, 

tj _ i-ct] ii) = hyvtlii) - h ivtj - \cf\ (0) = - h ivt\ - \ct] (0) , 

here we have (and we only need) central limit theorem for the sum of \ \yt\ — 
\Cf\ \ number of iid. a;i(0) variables, which proves the theorem. Similar argu¬ 
ment works for H < 0 also. □ 


5. The motion of the defect tracer 

With the help of another type of coupling, with any n G Z+, we prove L"- 
convergence for Q{t)/t of BL and totally asymmetric ZR models in this section. 
This coupling only works under convexity condition 1.5, which we assume for 
the rest of the paper. The idea of the proof is the following: we fix our (w, Q) 
pair and compare it with another model f. The difference between w and C is 
realized by second class particles. The current of these particles satisfies law of 
large numbers by separate ergodicity of w and and we compare their motion 
to our defect tracer Q placed on w. The main difficulty is finding the way 
to couple the defect tracer to the second class particles. As shown later, this 
coupling can not be made directly; we need to introduce a new process called 
the S-particles, a random process defined in terms of the second class particles. 

We set 0i < 02, then there exists a two dimensional measure /r on ZxZ, which 
has marginals and respectively, and for which p,{x, y) = O'li x > y. We 
fix two configurations ry and C of our model, distributed initially according to a 
product measure with marginals P(r]i{0) = x, 0 ( 0 ) = y) = p{x, y). Therefore, 
rj is itself in distribution , C, is in distribution , and 77 ^( 0 ) < 0 ( 0 ) for each 
site i is satisfied. According to the basic coupling described in subsection 3.1, 
Viit) 0 Ci{t) holds for all later times t, and we have a positive density of second 
class particles between these two models. The number of these particles at site 
i is fi — rji > 0. Hence they are initially distributed according to a product 
measure but, at later times, only the marginal distributions of 77 or of C will 
possess a product structure. Note that the joint distribution of the processes is 
translation invariant. 
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5.1 The Palm distribution 

For further applications, we want to select “a typical second class particle”. We 
do it as follows. We introduce the drifted form of the models: for fc € Z, 

(rfc k C)^ ■ — C^+fc’ 

li N € Z"*" is large enough, we choose uniformly one second class particle among 
the particles present at sites —N < i < N. We determine the distribution of the 
values of a function g depending on ( 77 , C), as seen from the position k of the 
randomly selected second class particle. For N large enough, the total number 

N 

j=-N 

of second class particles at sites —N < i < N is positive, and then 


{airk'n-, TkO) = E [E {girkV, TfcC) |^. C)] = 


N 


= E ^ giT.'q, 


Q - Vi 


Et-NiQ-Vj)J 

7+lJ2Z-N9{T^r], T,() ■ (Ci 


\i^-N 


= E 


2iV+l 


2N+1 J2j=-NiCj Vj) 



For bounded g, the random variable we see in the last line of the display is 
bounded, and is the quotient of two random variables, both having a.s. limit 
as ^ 00 by translation invariance and ergodicity of translations. Hence our 
expression converges due to dominated convergence, and have the limit 

(37) E {g{g, 0) : = lim E^ {g{Tk V, Tfe 0) = \ ' 

- N^oo - E(Cn-r?n) 


The distribution g defined by (37) is called the Palm distribution of the process. 
The Palm measure can be extended to non-negative functions g, see [17]. Note 
that P(Co (0) — 770 ( 0 ) > 0) = 1 according to this measure, i.e. we necessarily have 
at least one second class particle at the origin, if looking the process “as seen 
from a typical second class particle”. 

By initial product distribution of ( 77 , C)) A* is initially also a product measure, 
consisting of the original marginals fx for sites 7 7 ^ 0 , and of marginal 


(38) 


g{x, y) : 


g(x, y)- {y- x) 
E(Co - Vo) 


for site 7 = 0. For later use, we introduce the pair {ri'{t), Cit)) started from this 
initial product distribution /i. 


5.2 Random walk on the second class particles 

We label the second class particles between 77 and ( in space-order. Let (t) 
denote the position of the rri-th second class particle at time t. Initially, we look 
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for the first site possessing second class particle on the right side of the origin. 
We choose one of the particles at this site, giving it label m = 0: 

[/(o) ( 0 ) : = min {i > 0 : Q > r]i} . 

We label the particles at t = 0 in such a way that (0) < (0) (Vm G Z) 

(the order of particles at the same site is not important). We define J- ’{t) to 
be the algebraic number of second class particles passing the column between i 
and t +1 in the time interval [0, t\. This quantity is determined by the evolution 
of the processes rj and C,. For t = 0, we define 

(39) mi{Q) : = max{m : f7^'"^(0) < i\, 


while for t > 0 , 

mi(t) : = mi( 0 ) - \t). 

We label the particles at later times such that (39) holds at any time t as well. 
This method assures for all time t. The particles labeled 

from nii-i + 1 up to mi, exactly C,i — rji = mi — mi-i of them are at site i. (At 
sites i for which mi = mi-i, there is no second class particle). 

We have defined so far the coupled pair ry and C with the process of 

the second class particles indexed in space order at any time t. The latter will 
serve us as a background environment for a new random process, (s^”^(t))^gg. 
Initially, we put (0) : = n for each n. Assume that just before a second 
class particle jumps from a site t at a time t, £ {mi-i{t) + 1, mi-i{t) + 

2,..., TOi(t)}, which means '^(t) = i just before the jump. Then by the 
time t+ 0 of this jump, s("^(t + 0 ) : = ni(s(")(t)), where 11 ^ is a random uniform 
permutation on the integer set {mi-i{t) + 1, mi-i{t) + 2 ,... ,mi{t)}. 

We can represent this new process as follows. Initially, we put an extra 
particle, which we call S'-particle, on each second class particle. The S'-particles 
are labeled by n, and initially we put the n-th S'-particle on the n-th second 
particle. s^”^(t) stands for the index of the second class particle carrying the 
n-th S-particle. Whenever a jump of second class particle happens from site i, 
we permute uniformly and randomly the S-particles present at site i just before 
the jump. According to the labeling of second class particles, one jumping to 
the right (or to the left, respectively) from site i has index mi{t) (or mi-i{t) + l, 
respectively) and is carrying exactly the n-th S-particle, for which (^ + 0) = 
ni(s^"^(t)) = mi{t) (or mi-i{t) + 1, respectively). Hence a uniformly and 
randomly chosen S-particle is taken from the site i with the jumping second 
class particle. 

For simplicity, we define s{t) : = and S{t) : = [/('*(*)), and by simply 

saying the S-particle, we mean the zeroth S-particle at site S{t). Then S{t) 
represents a random walk moving always together with a second class particle, 
but having always probability l/(mj — mi-i) = 1/{Q — iji) of jumping together 
with a second class particle jumping from the site i. As can be derived from 
table 1 , the rate for a second class particle to jump to the left (or to the right) 
is f{—r]i) — /(—Ci) (or /(Ci) — /(^i)) respectively). Hence the rate for the S- 
particle to jump to the left (or to the right) together with the jumping second 
class particle from site i = S{t) is 


(40) 


f{-Vi) - /(-Ci) 
Ci - Vi 


(or respectively). 

Ci - v^ 
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Recall that S'(O) = (0) = is the hrst site on the right-hand 

side of the origin initially with second class particles. We introduce the notation 
W'it), Cit)) := (rs(o)^(t), Ts{o)C{t)), which is the C(0) process shifted 
to this initial position 5(0) of the 5-particle. We also introduce its 5"-particle: 
5"(t) : = 5(t) — 5(0). Hence the initial distribution of (?7"(0), C^^(O)) is modihed 
according to this random shifting-procedure; we show the details in the proof 
of the next lemma. 

Using the Palm measures, we show that the expected rates for 5 to jump 
are bounded in time. 

Lemma 5.1. Let n € k G Z, and 

(41) c,{t) : = f{Ci{t)) - f{r]iit)) + f{-Viit)) - 

the rate for any second class particle to jump from site i. Then 
E ([cs(t)(^)]" • [Csit)(t) - vs{t){t)f) < K{n, k) 
uniformly in time. 

Proof. First we consider the pair (??'(0), C^(0)) defined following (38). As de¬ 
scribed there, this is in fact the pair (r?(0), C(0)) at time t = 0, as seen from 
“a typical second class particle”, or equivalently, as seen from “a typical 5- 
particle”. In this pair, we have at least one second class particle at the origin, 
which we call S'. We let our process {rj', f') evolve, and we follow this “typical” 
5'-particle. Started from the Palm-distribution, this tagged 5'-particle keeps 
on “being typical” (see [17]), i.e. for a function g of the process as seen by 5', 

E {g{Ts^^t)r/{t), rs'(t)^{t))) = E {gig{t), C(<))) 
with definition (37). 

Now we first show the desired result for the 5'-particle of (ry', f') instead of 
the 5-particle of (ry, f). In the previous display, we put the function 

givit), Cit)) ■ = [co(i)]” • Koit) - Voit)f, 

and we denote by k~^ the positive part of k. We know that fo{t) — go{t) > 1 
holds P-a.s., hence 


E ([cs'(t)(0]” • - 7s'(t)(0]'') = E ([co(<)]” • [Coft) - voit)]'") < 

E([co(t)]"-[Co(0-7o(t)f^+^) 


< E 


([co(i)]”-[Co(t)-7o(t)]‘ 


E(Co(i)-ryo(t)) 


by (37). The function co(t) consists of sums of /(±ryo(t)) and f{±fo{t)), hence 
the numerator is an n -I- fc’*’ -I- 1-order polinom of these functions and of Co (t ), 
?yo(t). These are all random variables with all moments finite. Therefore, using 
Cauchy’s inequality, the numerator can be bounded from above by products of 
moments of either /(ryo(t)) or /(Co(0) or ?yo(t), or Co(0- The models g and C 
are both separately in their stationary distributions, hence these bounds are 
constants in time. The denominator is a positive number due to O 2 > 0i and 
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strict monotonicity of Eg( 2 :) in 6. We see that we found a bound, uniform in 
time for the function g of {rj', (^') as seen from S'. 

We need to find similar bound for a function g of the original pair ( 77 , ^), 
as seen from S. This is equivalent to hnding a bound for g of {rj", ({") dehned 
above, as seen from S" of this pair. Let us consider hrst the initial distribution 
of {rj", ("), which we shall call g". By definition, it is clear that this distribution 
is the product of the original marginals g for sites i > 0. Fix a K positive integer 
and two vectors x, y G Z^. For simplicity we introduce the notations 


%, 6] • = ■ ■ ■ • and j,] : = (Ca, ■ • •, Cb), 

X[a, b]-.= {xa, ■■■, Xb) and : = {ya, ■ ■ ■, Vb) 

and, where not written, we consider our models at time zero. We break the 
events according to the initial position 5'(0) of the S'-particle in the original pair 

iv, C): 


^ ilL[-K,0] -[-^. 0 ]’ i[-7f,0] y[-K,0]} 

= ^ illlSiO)-K, S{0)] =^[-i^,0]> ^![S(0)-/f,S(0)] ^^[-K,0]) ^ 

K 

= E P n] = 3i[-K, 0 ]. „] = 0 ]: 5(0) = n) + 

n=0 

00 

+ E P {iLln-K, n] = 0] > „] = 0] ’ ^(0) = ^) = 

K 

= EP t[n-iC.n] =^[-i^,0]. i[n-KM= 

n—0 

oo 

E P {ll[n-K,n] = ^l-K,0], C[r.-K,n] = ^[-iC.O]) ' fe, V) ' P{E„_k}, 

n—K-^1 

where the function En of x and y is an indicator dehned by 

En{x, y) ■■=1 {x-n = y-n, x-n +1 = y-n+ 1 , ■■■, X -1 = ?/_i, a;o < yo} , 
and the event En-K is 


Fn-K ■ = {VO = Co, »7l = Cl) ■ • ■ ) Vn-K-l = Cn-K-l} ■ 

The last equality follows from the product structure of g and from the fact that 
S'(O) is the hrst site to the right of the origin where rji yf Q. Continuing the 
computation results in 


P(!Z[-iC,0] -[-^.0]’ i[-7f,0] ^[-iC,0]) 


= n ■ 


i^-K 


K 


E ^"(^> y) + ^k{x, y) ■ ^ y{vo= Co}' 


I-K 


Ln=0 


n—K-\-l 


0 

= n ^(^*’ 2 /*) ■ 

i=-K 


■ K 

E y) + ^k{x, y) ■ 

.n—0 


y{vo = Col 
y{vo < Col 
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using translation-invariance. 

For later purposes, we are interested in the Radon-Nikodym derivative of 
the distribution /x" of ( 77 ", (") w.r.t. the Palm distribution of ( 77 ', Since 
both have product of marginals /x for sites z > 0 , we only have to deal with the 
left part of the origin. Passing to the limit K 00 , we have 


d^" 

dp 


(x, y) 


^(!Z[-Xf,o] -[-^.0]> C[-Xf,0] y.[-Kp]) 


— - - - 


■ K 

^ En{x, y) + Ek{x, y) ■ 

.n—0 


y{m = Col 
y{vo < Co} 


fJ'ixQ, yo) 

K^o, yo) 


00 

Enix, y) + lim Ek{x, y) ■ 
^ ~ K—*oc — 

.71 — 0 


y-ivo = Col 
y{m < Col 


yjxQ, yo) 
M^o, yo) 


■ y) 

71—0 


for /{-almost all configurations {x, y). Note that the sum on the right-hand side 
gives exactly the distance between the origin and the first position i to the left of 
the origin with Xi ^ yi- Hence this sum is finite for /t-almost all configurations 

fe y)- 

In view of this result, we can now obtain our estimates. The main idea here 
is that the pairs ( 77 ', CO iv"^ C") only differ in their initial distribution, 

hence their behavior conditioned on the same initial configuration agree. This 
is used for obtaining the third expression, and Cauchy’s inequality is used for 


31 



the fourth one below. 


E {[cs"(t)it)r • [Cs"{t){t) - ils"(t){t)]'') — 

= y E ([cs„(i)(t)]" • [Cs"{t){t) - \r/' (0) = X, C{0) = y) X 


nn{a:o<yo} 


X y) = 

E([cs/(t)(t)]” • [Cs'{t)it) - vs'it){t)f \rfiO) =X, C'(0) =y) X 


On{a:o<yo} 


X ^ En{x, y) d^(x, y) < 

y{xo,yo) - - - 

I [E ([cs,(,)(t)]" • [Cs'(,)(t) - vs'it)it)]’^\n:i0)=x, C'(0) = y)]' X 

slnticoXiyol 


X dy{x, y) 


n{a:o<yo} 


y{xo, yo) 

M^o, yo) 


^ En{x, y) ■ dg(x, y) < 


I E i[cs>it){t)r ■ IMt) - ys'it)it)r\v:iO) = X, C'(0) = y) X 
T2n{a:o<i/o} 


X dyix, y) 


yjxp, yo) 
Kxo, yo) 


n{a:o<yo} 


X! En(x, y) ■ d^(x, y) = 


= [E {[cs-^t){t)r ■ Ks'wit) - vs'it)it)r)] " X 


lr2n{xo<yo} 


E(Co - Vo) j:, r ^ j ^ ^ 

-> En{x, y) ■ dfx{x, y) 

yo - ^0 L^o “ J “ “ 


by (38). The first factor of the last display is finite by the first part of the proof. 
Using the definition of the indicator E„, the second factor can be bounded from 
above by 

1 

2 

/ r oo “I 

^(2n + 1) • En{x, y) ■ d^(x, y) = 

Ln^O J 

nn{afo<yo} 


LXJ 

= [E(Co - Vo)]^ ■ ^(2n + 1) • ii{r]o = Co}” • y{vo < Co} 


using the product property of yt, and is again finite. 


Using the rates for the S'-particle to move, we can prove the following bound 
for the moments of S{t): 
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Proposition 5.2. For n € Z+, 


(42) 



< K(n) < oo 


for all large t. 

Proof. For this proof, we denote the jumping rates (40) for the S'-particle by 
^Sieft right ^ respectively. For t > 0, we consider the derivative of the 

quantity above, using these rates: 

= _^E(|S(()r) + l^E(|S(()l") = -^E (15(01") + 
+ 1e • (|s'(t) - ir - is'(t)r)+r^"s“(-i5'(t) + ir - is'(t)r)i. 

t” L J 

For |S'(t)| > 1, we can bound our expressions: 



< - 


t"+i 


E (1^(01”) + —E ((r^"s“ + • |^(t)r-^) . 


We continue by using Holder’s inequality on the right-hand side: 


(43) 


dt \ t" 


< 


< -E 

t 


V t” 


^ {E [C 


.S right 


left\n 




E 


V t" 



Recall that 


right (^) ^ Cs(t)(t) • [Cs(t) (0 - ?7S(t) (<)] \ 

hence lemma 5.1 is applicable with k = —n to show that 

E [(r®"s*'*(t) -hr^^®f‘(t))"] 

is bounded in time. Therefore, (43) can be written in the form 



< -E 

“ t 



K'{n) 
+ — 



with some positive constant K'{n). This means that E(|5'(t)|"/t”) is bounded 
from above by a solution of the differential equation 


y{i) 


n , , K'(n) , , 1 


Observe that the right-hand side is negative whenever 


yit) > 



hence assuming y{to) < oo for some to > 0, y{t) is bounded (for all t > to), 
which gives the proof. □ 
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Now we show law of large numbers for s{t), and then we can show law 
of large numbers for S{t). For what follows, E' stands for the expectation 
values according to the distribution of ry, i-e- our background 

process which determine also. Let J^{t) denote the cr-field containing all 

information about these quantities at time t. Then J^{t) contains all randomness 
except for the random permutations on With (41), we also introduce 

the notations 


C'i(t) : = {mi(t) - ■ Ci{t), 

p{y, t) : = P(s(t) =y \ T{t)), 

(44) 

and 


= max p{y,t)- min ^p{y,t) 

Tni — i{t)<y<mi{t) mi-i[t)<.y<Tni[t) 

if mi{t) — > 1, and Ai{t) : = 0 otherwise. 

Lemma 5.3. 


(45) 


dt 


E(|s(i)l) < 


\ 


E' ^ A{t) 


\ 


E' i; 


Proof. We use convention that the empty sum equals zero. 


,iE(|,(,)|) = lim E(Wt + .)l)-E(|.W|) 

dt e^o e 


= lim > 
e^O 

Z — — 00 

oo mi{t) 


P(s(t + e)=z)- |z| - P(s(t) = z) ■ \z\ 


P(s(t + e)=z \ T{t)) ■ |z| - P(s(t) = z I T{t)) ■ |z| 


i= —oo z=mi-i(t) + l 


We know that uniform random permutation on the indices present at site S 
happens at each jump of second class particles from i at time t. The basic idea 
is that this permutation makes the probabilities p{y, t) = P(s(t) = y\ T{t)) 
equalized between y = + 1... TOi(t). This jump happens with rate Ci{t) 

defined in (41), hence for a site i with at least one second class particle and for 

TOi_i(t) + 1 < z < 


P(s(t + e) = z I T{t)) 


{1 - eci{t))P{s{t) =z| J^{t))+ 


rrii (t) 

+ eci{t) ^ 


Pis{t)=y\ Ht)) 
mi{t) - 


0(e). 


Then we obtain 


dt 


E(|s(t)|)=E' 


E 


rrii (t) 

.w i: 

z—mi — i (i) + l 


rrii {i) 

E 


p {y, t) 


mi{t) - TOj_i(t) 


-p{z, t) • 
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There exists a tt^ permutation of the numbers {mi-i{t) + 1... mi{t)}, for which 


mi{t) 

E E 


piy, t) 


unit) - 


TTLi {t) 

z\< ^ p{z,t)-\7ri{z)\ 

z—rrii-i (t) + l 


holds (by permuting higher values of \z\ on higher weights), and hence 


d 

dt" 


rrii (t) 


-E(|s(t)|) < E' ^ Ci(t) p{z,t)-i\Tr,{z)\-\z\) = 

i— — oo z—mi-i{t)-\-l 

oo mi{t) 

= E'Ec*(0 Y B(y, t)] • (|7r,(z)| - |z|) < 


i=-oc z=mi-i{t) + l 


i(t) 


<E' ^ c*(t) Y ( p{y,t)- 

- min p{y,t)] ■ - mi-i{t)) = 

= E' V c*(t)( max p{y,t)- 

\mi_i (t)<v<mi(t) 


i : rrii (t) >mi_i (t) + l 


oo 

- min p{y, t)) ■ {mi{t) - =E' Y Ai{t)C^(t) 

2 = —CX3 

with definitions (44). Finally, we use Schwartz and Cauchy’s inequality (for 
simplicity we do not denote time-dependence of the quantities below): 


dt 


E(|s(t)|)<E' Y A^Ci=^' Y 


< E' 




E 


J^-OC 


< 




E' E 


j = -oo 


□ 


Lemma 5.4. The expression 


\ 


E' E A,(t)Cf(t), 


j = -co 


which is the second factor on the right-hand side of (45), is a bounded function 
of time. 

Proof. Due to definitions (44), Ai can be bounded from above by 

rrii (i) 

Ai{t)< max p{y,t)< Y Piy,t) = PiSit) =i\ :F{t)), 

y—mi-i{t)-\-l 
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the probability that our S'-particle is at site i. Hence 


\ 


E' ^ 


\ 


E' ^ p{sit)=j\m)-c]{t) = 


1 


\W 

w 




The expectation in the last term is bounded in time by lemma 5.1 with n = 
2, k = 4, since 

C\ = {m^ - mi-if' • Cj = Cj • (Cj - rjif . 


□ 


As we know, for any site i, the probabilities p(j/, t) can only change by 
equalizing between y = mi-i{t) + and the initial distribution is 

concentrated on {s{t = 0) = 0}. Therefore, at every moment t, the function 
y p{y, t) is unimodal. This is clearly the initial situation, and it stays true 
after each change of this function. By the equalizing property of the (p (y, 
process at a jump of second class particle from site i, 


max p (z, t) 

mi-1 {t)<.z<mi (t) 

can never increase. Hence the global maximum max 2 gzp(z, t) is also a non¬ 
increasing function of t, and it is bounded as well. Thus its limit exists, which 
we denote by P. It is believed that P = 0 but we cannot prove this, and this is 
not necessary for our arguments. 

Lemma 5.5. Assume P > 0. Then the set 


{x G Z : p{x, t) > P} 

is always contained in the interval [—1/P, 1/P]. 

Proof. The statement clearly holds initially. For a discrete interval [x, y] (with 
possibly X = y as well), we introduce the block-average 

1 ^ 

y] (^) ■ =-TT P 

^ ^ y — X + I 

and we say that [x, y] is a good block, if < min^gfa, j,] l/\z\ (for site 

z = 0, we can write 1 instead of l/|z|). Any interval is a good block initially. 
We show this for any time t as well. More precisely, fix a; < y, and assume that 
at a moment t, an equalization in the interval [rt, ri] happens: 

p{z, t + 0)= P[u,t,](0 

for each 2 ; G [rt, ?;]. If each finite interval is a good block at t, then we show that 
[x, y] is also a good block after this step, at t -I- 0. There are four cases. 

(i) If [u, u] and [x, y] are disjoint or [u, v] C [x, y], then the block-average of 
[x, y] does not change by this step, hence it keeps on being a good block. 
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(ii) If [x, y] C [u, v], then + 0) = + 0) = and [u, v] 

was a good block at time t, hence [a;, y] is also a good block after this step. 

(iii) In case [x, y] \ [u, v] ^ 0, [u, v] \ [x, y] ^ ^ and 
before the step, then 

-^[•u, f] “b 0) = ii] (t) ^ B^x,y]\[u,v]i^^ ~ B^x, y]\['U, i?] (t “b 0), 

hence B[x,y]{t + 0) < B[x,y]ulu,v]{t + 0)- The latter does not change by 
the step, thus [x, y] U [u, w] keeps on being a good block, which shows that 
[x, y] is also a good block after the step. 

(iv) In case [x, y] \ [u, v] ^ 0, [u, v] \ [x, y] 7 b 0 and B[„,„](t) < y]\[«,«] (0 

before the step, then by unimodality, i?[„ „](t) < since the 

function z —>■ p{z, t) has no local minimum. This means that t)]n[a:, y] 
does not increase: 


B[u^v]n[x,y]{^ “b 0 ) “b 0 ) ^ i!]n[a:, y] (^) • 

Since B[x,y]\lu,v] does not change, Bix,y] can not increase either, and [x, y] 
was a good block before the step, thus it keeps on being a good block. 

Applying this result shows the interval containing any single point z to be a good 
block, i.e. p (z, t) <1/P for z ^ [—1/P, 1/P], which completes the proof. □ 


Lemma 5.6. Assume limj^oo (z, t) = P > 0. Then there are z, y 

neighboring sites in the interval [—1/P — 1, 1/P + 1] and a time T > 0, such 
that the second class particles indexed by z and y cannot he at the same site 
after T: U'-^\t) 7 b [/(y)(t) (Vt > T). 


Proof. Let 


A : = < z G : lim supp ( 2 , t) = P 


^ 0 - 


By the previous lemma, A C [—1/P, 1/P], and any index Zmax(t), for which 
P (zmaxit), t) is maximal (and hence larger than or equals to P), is also contained 
in [—1/P, 1/P] for any t. With fixed Pi < P large enough, there exists a 
moment T, such that p (x, t) < Pi for any x ^ A and for all t > T. Hence 
by p(zmax(t), t) > P, all indices Zmax(t) G A for all t > T. Let us fix z G A 
and y ^ A neighbors, and y' (/ A the other neighbor of A. Then infinitely 
often for t > T, p {z, t) > P > Pi > p{y, t) and Pi > p{y', t) happens. 
In this situation, assume that p (z, t) decreases due to equalization with its 
neighbors in A. Would the result of this step be p (z, t) < P, all indices Zmax(t) 
would be included in this step by unimodality, hence p(zmax(t + 0 ), t + 0 ) < 
P would follow, a contradiction. Thus we see that p (z, t) > P can only be 
violated by an equalization including y or y'. If this equalization also includes 
all indices Zmax(t), then the result must be p{y, t) > P or p{y', t) > P hy 
P {zmax{t + 0), t + 0) > P, pulling out at least P — Pi probability from the set 
A. If this step does not include all Zmax indices, then it includes indices all 
with probability at least P, hence pulling out at least (P — Pi)/2 probability 
from the set A. Since t > T, Zmax(t) G A, and hence by unimodality, the 
joint probability of the set A can only decrease. We conclude that assuming 
equalizing of probabilities between z G A and y or y' (/ A infinitely often results 
in decreasing the joint probability of the finite set A infinitely often by a positive 
constant, which contradicts P < p (zmax(t), t) and Zmax(t) G A. □ 
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Now we can prove law of large numbers for the index s{t) of the second class 
particles carrying our S'-particle: 

Proposition 5.7. 

(VJ > 0) ^Ihn P ^ >(5^ = 0. 

Proof. By the previous lemma, we see that for P > 0 there exists a neighboring 
pair {z, y) G [—1/F — 1, 1/P + 1] of second class particles which will never 
meet after some T. After T, the process s{t) can not cross such a pair (z, y). 
By translation invariance, it follows a.s. that such pairs appear with positive 
density on Z in this case, thus s{f) is bounded a.s. and the statement is true. 
Hence we assume P = 0 for the rest of the proof. By unimodality. 


OO 

(46) ^ A,{t)= ^ ( max p{y,t)- 

i— — oo i : Tni{t)'>7ni — i{t)-\-l 

— min p{y,t)) <2 maxp( 2 ;, t). 

1 (t)<y<mi (t) / zGZ 


Indirectly let’s assume 


(35 > 0) {3K > 0) (VT > 0) {3t > T) : P 


>6] > K. 


Then it follows that 


(47) E(|s(t)|) > 

for infinitely many and arbitrarily large t > 0. By (46) and P = 0, 

OO 

i— — oo 

thus by dominated convergence theorem 

0 . 

Hence by lemma 5.3 

, OO OO 

-E(|s(t)|) < . E' E Mt) • . E' E ^ 0 

\ i— — oo \ j — — oo 

when t —>■ OO, as 

OO 

E' E 

\ j^-oo 

is bounded by lemma 5.4. That means that 

|e(|.(,)|) 

tends to zero as t ^ oo, which contradicts (47). □ 
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Now we show the law of large numbers for S(t), the random walk on the 
background process 77 with parameter 61 and C with parameter 62 . 

Proposition 5.8. Let 


(48) 


c{ 0 i, 62 )-.= 2 


cosh( 02 ) — cosh(0i) 

E02(Co) - Eej(?7o) 


for BL models, and 


(49) 


c(^i, 6 * 2 ) 


Ee2(Co) - ^eAvo) 


for the ZR process. Then for every (5 > 0 


(50) 


lim P 

t—*oo 


t 


c(6*i, 6*2) 



= 0 . 


Proof. We show the proposition for BL models, the modification for the ZR 
process is straightforward. By the coupling rules, if a second class particle 
jumps from t to z + 1 then the column gi of C between sites i and z + 1 increases 
by one. If one jumps from z + 1 to z then the column hi of 77 increases by 

one. Hence for the current Jl of second class particles defined earlier in this 
subsection, 

JT \t) = - ffi(O)) - {ht{t) - /i*(0)), 

i.e. it is the difference between the growth of columns z of C and of 77 until time 
t. Due to separate ergodicity of each f and 77 , we have law of large numbers 
for gi{f) — gi{0) and for hi{t) — hi{0), since each of these models is distributed 
according to its ergodic stationary measure. Hence with the expectation of the 
column growth rates, we have 


(51) 


lim 

t^OO 



t 


E«2 (/(Co) + /(-Co)) - E,, (/(770) + f{-vo)) = 

= 2 (cosh( 6 l 2 ) — cosh( 0 i)) a.s. 


We extend definition (39) for a; G R: 

iTT-xit) := m\x\ (t) = max{? 7 z : U^'^\t) < x} 

Obviously, m^it) = mx{0) — J‘^\^{t). If RT G M then 

lim = E(Co) - E( 77 o) = : p a.s. 

v^oo /I y 

since at t = 0 , the starting distribution of the number of second class particles 
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at different sites is a product measure. 


(52) lim P 

t—*oo 


= lim P 

t —^OO 


= lim P 

t —^OO 


^iKt){t) ^ e2)-K)p >e^ = 

p ^ ^ — Kp + 2 cosli(d 2 ) — 2 cosh(0i) > = 


= lini p 


— ifp + 2 cos 1 i( 02 ) — 2cosh(di) >e = 


+ 2 cosli( 02 ) — 2 cosh(di) > e = 


= lim P — 

t—^OO \ 


—2 —^—- + 2 cosh(d 2 ) — 2cosli(di) >£1=0 


by translation-invariance and by (51), for any £ > 0. Recall that S{t) is the 
position of the zeroth S'-particle, i.e. the position of the s(t)-th second class 
particle: S{t) = Hence 


(53) P ( ^ - c(0i, 02) > ^) = P ( ^ = 

(u^sit)) , 02 ) t + <51) + p (t) < c(0i , 02 ) t-St 


In case 

[/(®(‘»(t) > C(01, 02)t + '5f 

it follows by definitions of (t) and of p that 
s{t) > m(^ci 9 u 92 )t+st){t) , hence 

P ([/(®(*»(t) > c(0i, 02 )t + St^ < P 


<p(f >^d+p 


'>n(c(9i,92)t+5 t) (t) ^ A 

t ^ 2^/ ■ 


As time goes on, the first term goes to zero due to proposition 5.7, and so does 
the second term by (52) (with K = c(0i, 62 ) + S). 

In case 

[/("(‘»(t) <c(0i, 02 ) 1-81 

it follows that 

s{t) < m(ci 9 u 92 )t-st){t) , hence 

P (^{7(^(*»(t) < c(0i, 02 )t- 5 t^ < P < 


p(i^<-|p)+p 


m(c{0i,92)t-st)it) 8 

—^^— > —ttP 
t 2 


The first term again goes to zero due to proposition 5.7, and so does the second 
term by (52) (with K = c(0i, 02) — 8 ). Thus we see that both terms on the 
right-hand side of (53) tend to zero as t ^ 00 . □ 
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5.3 Coupling the defect tracer to the ^-particles 

We fix the model w in stationary distribution with the defect tracer Q{t) 
started from the origin. We prove theorem 1.6 for BL and (totally asymmetric) 
ZR models. A natural idea would be to couple the defect tracer Q to the second 
class particles, present at the same site Q. The problem is that, either to the 
left or to the right, the rate for any jump of second class particles form the site 
Q may be higher than the rate for Q to jump. On the other hand, one second 
class particle always stays at site i after one jump from i, in case more than 
one of them were present at i. The solution is to couple the defect tracer to the 
S'-particle, for which the desired conditions are already proven by propositions 
5.2 and 5.8. For simplicity reasons, in case of the ZR process we let f{—z) : = 0 
for 2 ; > 0, and hence of ZR is also zero in these cases. 

The upper bound for Q. 

First, we identify ij distributed according to /i^ with w possessing the defect 
tracer Q, therefore we set 6*1 := 9 < 02 - We have then uJi{t) < (i{t) for all 
t according to the basic coupling, and recall that Q(0) = 0 < S'(O). In what 
follows, we are going to couple the random permutations of the S'-particles, thus 
the random walk S{t) of the zeroth A-particle, with the defect tracer Q{t). We 
only couple them in case Q{t) = S{t). The basic observation we use is that 
the rates (40) for the jump of the S'-particle can be compared to the rates for 
the jump of the defect tracer Q{t). As we have seen at the introduction of 
BL models, it is enough to consider the “effect of bricklayers” standing at each 
position i. That is to say, we are allowed to consider the w^-dependent parts of 
r{u)i-i, tOi) and r(a;j, w^+i) only, since the Wi-dependent parts are added to the 
Wi-i-dependent or to the w^+i-dependent parts in these rates. In the rest of the 
paper, we describe couplings by giving rates of bricklayers standing at each site 
i. This observation also holds for the zero range process (by saying rate for a 
particle to jump instead of saying rate for bricklayers to lay bricks). 

In tables 2 and 3, hi ] means that the column of the model w between i and 
i \ has increased by one, gi ] means that this column of C has increased by 
one, r\ means the jump to the right from i, means the jump to the left from 
i. 

Note that by z = S' = Q, Q > oJi + 1. The rates are non negative due to 
monotonicity of / and convexity condition 1.5. By summing the rates corre¬ 
sponding to any column of the tables, one can verify that each w and C evolves 
according its original rates, Q has the jump rates according to the basic coupling 
described in table 1, and S also has the appropriate rates (40). We see that 
once being at position S, the defect tracer Q can’t move right without moving 
S with it and S can’t move left without moving Q with it. Hence our rules 
preserve the condition Q < S. 

We have so far the upper bound Q{t) < S{t), and we have the law of large 
numbers (50) with speed c{ 6 , 62 ) defined in either (48) or in (49) for any 62 > 9, 
and the n-th moment condition (42) for this S{t) process. 

The lower bound for Q. 

Now we show a similar coupling which results in a lower bound for Q. The 
natural idea would be to identify C with w, and couple Q to the S'-particle. The 
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rates for Q and S to jump with would allow Q{t) > S{t). However, this coupling 
can not be realized in a similar way that the coupling described above: there is 
no way for Q and S to step together, since only one brick can be laid at a time 
to a column. 

Therefore, we need to modify the initial distribution of the models as follows. 
Let ^{x, y) be, as before, a two dimensional distribution giving probability zero 
to X > y, and having marginals /igj and respectively. Fix the pair {rj, C), 
as before, with the product of fj,{x, y) for different sites as initial distribution. 
Define 

(54) y!{y, x) : =/x(x, y) ■ 

Fix the pair ( 77 ', C^), with the product of /i(a;, y) for each site i 7 ^ 0 and of 
y'{x, y) for the site i = 0 as initial distribution. Then 77 '(0) < Ci(0) holds a.s. 
for each site i, hence the basic coupling is applicable for this pair of models. We 
have second class particles between 77 ' and and we introduce the S"-particles 
as well, starting S'g from the first site on the left-hand side of the origin: 

S'{Q) = 5'(0) := max {7 < 0 : C'(0) > r;'( 0 )}. 

Assume now that the S'-particle of ( 77 , C) is also started from the first site on 
the left-hand side of the origin, instead of starting it from the right-hand side 
of the origin: 

S'(O) = S'o(O) : = max{i < 0 : Ci(0) > ’7*(0)}- 

Then it is clear, that propositions 5.2 and 5.8 also hold for this S'-particle. Now 
we derive these statements for S' as well. Since initially ( 77 ', (^') only differs from 
( 77 , C) by the distribution at the origin, the conditional expectations 

(55) E (S'(t) 1 77 '(0) = x, Co(0) = 2/)) = E (S(t) 1 770 ( 0 ) = x, Co(0) = y)) 
agree. This is the basic idea of the following 

Lemma 5.9. The moment condition (42) and the law of large numbers (50) 
hold for S' as well. 

Proof. By the use of (55) and Cauchy’s inequality in a similar way than in the 
proof of lemma 5.1, 



1 

2 
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The first factor of the display is bounded by proposition 5.2. For the second 
factor, by (54) and ( 6 ) we write 




x<v 


x<y 


^^eAy) 


_ ^^2 {y ~ 1) 

(y) 




fiy) 


yGZ 


e^= f{y-1)1 Z{ 02 ) e 


^ =4-E.,(/(2/)^), 


which is again finite. Hence (42) holds for S' as well. 

For the law of large numbers, we know that for any 5 > 0, 


0 = lim P 

t^OO 


Sit) 


— c( 0 i, 6 * 2 ) 


>S] = 


= lim V P 

f. —^r>o I J 


x<y 


Sit) 


— c(6*i, 62 ) 


>6 


? 7 o( 0 ) = X, Co(0) = y] ■ yix, y) = 


= lim P 


x<y 


S'it) 


— c( 0 i, 6 * 2 ) 


> .5 


»7o(0) =x, Co(o) =J/) ■ yix, y), 


hence (50) follows for S' as well by absolute continuity of y' w.r.t. /x. 


□ 


In order to obtain lower bound for Q of w distributed according to /i^, set 
02 = 0 > 01- The marginal distribution of Co(0) is the second marginal of /x', 
namely, iy— 1) = yeiy— 1). Hence it is possible to fix the pair iy', C') defined 
above with 

Cit) = fti(i) + ^Q(t) j Q(o) = Oj 

i.e. w is coupled to C' with the defect tracer Q between them. Note that S''(0) < 
0 = <5(0). We show the coupling that preserves S'it) < Qit) for all later times. 
We only couple Q to the random permutations acting on S' in case Q = S' 
for a site i. For tables 4 and 5, /i' t means that the column of the model y' 
between i and x + 1 has increased by one, t means that this column of C,' has 
increased by one. Note that by z = 5" = Q, Cx > z/' + 1. As at the coupling 
for the upper bound, the rates are non negative due to monotonicity of / and 
convexity condition 1.5. By summing the rates corresponding to any column of 
the tables, one can verify that each y' and C,' evolves according its original rates, 
Q has the jump rates according to the basic coupling described in table 1 (hence 
w also evolves according its original rates), and S' also has the appropriate rates 
(40). We see that once being at position S', the defect tracer Q can’t move left 
without moving S' with it and S' can’t move right without moving Q with it. 
Hence our rules preserve the condition Q > S'. 


Proof of theorem 1.6. By the upper bound and the lower bound above, we have 

Sit) > Qit) > S'it) 

and for any 02 > 0 > 0i, have weak law of large numbers for S with c(0, 02), 
and for S' with c( 6 *i, 0 ), respectively. Hence taking the limits 0i Z' 0 and 02 \ 
completes the proof of the law of large numbers ( 10 ) by computing 

lim c( 0 i, 0) = lim c( 0 , ^ 2 ) = ^( 0 ) 
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both for BL and ZR models. Moreover, for any n S Z+, we have n-th moment 
condition (42) for both S and 5", hence not only (11), but the n-th moment 
condition follows as well for Q. This also shows ^"-convergence of Q{t)/t for 
any n G □ 

5.4 Strict monotonicity of C{9) 

As a consequence of the type of coupling methods shown above, we are able to 
show strict convexity of the function Ti.{g) of (8). First we refer to the coupling 
which shows (non strict) convexity, and then we complete the proof of strict 
convexity by some analytic arguments. 

Remark 5.10. Let gj_, u/ be two copies of a model (either BL or ZR model) 
possessing condition 1.5, with the defect tracers Q{t) and Q'{t), respectively. 
Assume that for each site i and for time t = 0 

Wi(0) < u;'(0) and <5(0) < <5^(0). 

Then it is possible to couple such way that for allt>0 and any i € I, 
uji{t) < uj'i{t) and Q{t) < Q'{t) a.s. 


is satisfied. 

This coupling is very similar to the ones shown in this subsection, we do not 
give the details here. The pair (w, u/) is coupled according to the basic coupling, 
and we can apply this proposition for the case when their joint distribution has 
marginals pL^ and , respectively. Then we simply see that the motion of the 
defect tracer of a model has a monotonicity in the parameter 9 of the model’s 
stationary distribution. In the introduction we saw that this implies convexity 
of the function Tt{g). We prove now strict convexity of this function: 

Proof of proposition 1.7. First note that by the form (6) of the measure pig, we 
have 

gi9)= E,(u;) =^log(Z(0)), 

Ee (^2) = ^Ee (a;) > 0, 

E, P) = ^E, P) = ^(Var,(o;)), 

where tilde stands for the centered variable. For the BL model, we need to show 
strict convexity of the function 

H(p)=E,(,)(r)=e''('?)+e-^(e). 


We compute its derivative 


de ^ ^ E,(S2) ’ 


and, similarly, the second derivative 

j2 


-hnit) = ^ 


dg^ 


[E, 


[(e^ + e-^) E, (S52)_(e«_e-^) E, . 
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Hence (strict) positivity of 

(56) [(e^ + e-^) E, _ (e^ _ e"^) E, 

on an interval of 6 is equivalent to (strict) convexity of 7i(p) on the corresponding 
interval of g{9). (56) contains derivatives of log (Z (0)), which is by definition 
analytic, hence (56) is also an analytic function of 6 . Moreover, by the previous 
remark, we know convexity of hence non-negativity of (56). Since this 

function is strictly positive at d = 0 by symmetry properties of /ig, there are 
at most countably many isolated points at which this analytic function is not 
strictly positive, hence we have at most countably many isolated points at which 
the second derivative of 'H{q) is not strictly positive. This completes the proof 
for the BL models. 

As for the ZR process, similar computation leads to 

[e^E, (152)-e^E, 

in place of (56). As we know non-negativity of this function by convexity of 
'H{q), we only need to show Eg (w^) ^ Eg (w^) for some d, then the previous 
analytic argument leads to strict convexity. 

Indirectly, assume 


(57) 


Eg (a;2^ = Eg (w^) 


for all 0 < 0. Since the right-hand side is the derivative of the left-hand side, it 
follows that 

E, (152) = A. e^ 

for some A > 0. Integrating this we have 

E0(w) = A-e^ 


(the additive constant is zero as can be seen by taking the limit 0 —> —oo). 
Integrating again we have 


log {Z{9)) = A • e® -b AT, i.e. 
Z(6») = AT'-e^®", i.e. 


OO 


E 



K' 


E- 

z=0 


• e 


ez 


z! 


for all 0 < 0, which leads to f{z)l = zlfA^, f{z) = zjA. Hence we see that if 
at least for one z > I value we have f{z + I) — f{z) > f{z) — f{z — 1), then 
(57) is not true for some 0, and then strict convexity of (p) holds. We also see 
linearity of 7i(p) when / is linear. □ 
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Q ^ 
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a second class particle 

fi-iVi - 1) - f{-Q) 

• 





[/(-‘^i) - /(-Wi - 1)]- 

Ci-^i 

B 

1 

H 

1 


Ci-“i 

• 


• 

• 


/(-G) 

• 

• 





Table 2: Rates for Q and S to step left and for bricklayers at site i = S = Q to 
lay brick on their left 
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with rate 
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Q r\ 

S rx 

a second class particle 

X [/(Ci) - 

1 

B 

1 
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/(Ci)-/K) 

- [/(Wi + 1) - 

1 

B 

1 

H 


f{uJi + 1) - f{uJi) 


• 

• 

• 

rv 

/(w*) 

B 

B 





Table 3: Rates for Q and S to step right and for bricklayers at site i = S = Q 
to lay brick on their right 
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C-iT 

C-iT 

Q ^ 

S' ^ 

a second class particle 

Ci-Vi 

■ 

B 

1 

1 



B 

B 

1 

B 


/(-C' + i)-/(-CD 

• 


• 

D 


/(-CO 

• 

D 





Table 4: Rates for Q and S' to step left and for bricklayers at site i = S' = Q 
to lay brick on their left 
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